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Abstract 

The Rabi Hamiltonian, describing the coupling of a two-level sys- 
tem to a single quantized boson mode, is studied in the Bargmann- 
Fock representation. The corresponding system of differential equa- 
tions is transformed into a canonical form in which all regular singu- 
larities between zero and infinity have been removed. The canonical or 
Birkhoff-transformed equations give rise to a two-dimensional eigen- 
value problem, involving the energy and a transformational parameter 
which affects the coupling strength. The known isolated exact solu- 
tions of the Rabi Hamiltonian are found to correspond to the uncou- 
pled form of the canonical system. 

PACS numbers: 03.65.Ge, 33.20.Wr. 



I. Introduction 



In the study of dynamical problems the harmonic oscillator occupies a promi- 
nent place, as a prototype of the fundamental unitary symmetry group. The 
spell of group theory also extends to anharmonic oscillators, which have re- 
cently been exposed as mere q-deformations of the unitary group algebra |I| . 
In contrast dynamical problems which arise in the study of boson-fermion 
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interactions are more reluctant to reveal their hidden symmetries. Such prob- 
lems are more involved since the actual states of the system are of composite 
nature with a boson and a fermion part. This results in the appearance of 
singularities in the corresponding eigenvalue problem. Examples include the 
Jahn- Teller Hamiltonian in molecular physics, and the Rabi Hamiltonian in 
nuclear physics and quantum optics. 

A rather peculiar feature of these systems is the emergence of isolated 
exact solutions, with eigenvalues that correspond to simple expressions in 
rational numbers. Such were obtained - in a rather heuristic way - by Judd 
|2| for the case of the E £g> e Jahn- Teller Hamiltonian, and by Kus || for 
the spectrum of a two-level atom coupled to a single quantized mode. It 
was suggested that the exact solutions probably hint at some dynamical 
symmetry group - but so far no progress was reported in this direction. 

In the present paper we analyze the one mode Rabi problem in a more 
rigorous way. First - following || and [|J - the dynamical problem is defined 
in the Bargmann-Fock Hilbert space of entire functions. The resulting system 
of differential equations is then put into canonical form using a theorem due 
to Birkhoff. Under this transformation the Kus exact solutions are found to 
be mapped onto the levels of a displaced harmonic oscillator. In this way 
hidden symmetry appears. 

II. The Bargmann-Fock Representation of a 
Hilbert Space 

Before turning to the actual results we briefly review basic information about 
the Hilbert space of entire functions introduced by Fock || and Bargmann 
0. Let J-'n, for n integral, be the set of entire analytic functions f(z), where 
z = (zi, ...,z n ) and Zp. 6 C are complex numbers, k = 1, ...,n. Because f(z) 
is entire it has an everywhere converging power series 

/0)= E asw^i 1 •••■•4™> (!) 
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where summation extends over the whole set of non-negative integ ers ki, . .., k n 
and OLk 1> ...,k n are complex coefficients. We define an inner product of two ele- 
ments / and g of T n by 



{f, 9) = / f( z ) ■ g{z)dHn{z) 



(2) 



where f(z) ■ g(z) is a usual scalar product in C", 



dfin(z) = — exp(-z ■ z) JJ dx k dy k , 



TV 



z k = x k + iyk 



(3) 



fe=i 



and the integration extends over the whole space C n . The Bargmann-Fock 
space jF n is a set of all entire functions (H), which have a finite norm (/, /) < 
00. This is equivalent to the requirement that 



\f{z)\ < cexp(-7z-z) 



(4) 



where c and 7 are positive constants with 7 < 1. The Bargmann-Fock space 
T n with the inner product defined by (0) is a Hilbert space. 

Let us now consider two operators in T n : multiplication by z^ and dif- 
ferentiation Since the functions f(z) of T n are analytic, Zfc/ and -^-f 
always exist. The operators satisfy the commutation rules 



[zk, zi] = 0, 



d d 

dz k ' dzi 



0. 



_d_ 

dzk 



zi 



Ski. 



(5) 



Furthermore, with respect to the inner product 
conjugate 

(Zkf,g) 



Zk and -r- are Hermitian 



> dz k , 



(6) 



whenever Zkf and belong to JF n . 

On the other hand relations @ and are well known algebraic relations 
defining Hermitian conjugate annihilation and creation operators of 
boson fields in second quantization, 



a k, a i 



0. 



[a k ,ai] = 0, 



>kl- 



(7) 
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We conclude therefore that within T n the annihilation operator is repre- 
sented by the operation and the creation operator a\ corresponds to the 
multiplication by z^. 

As an instructive example let us take a set of n identical uncoupled har- 
monic oscillators which are described up to a constant by the Hamiltonian 
H = J2?=i a t a i- The corresponding operator in the Bargmann-Fock space 
(denoted by H) is TC = Y^=\ zij^ and the corresponding eigenproblem is the 
system of differential equations 

n d 

E*J-fk(*) = E k fk(z). (8) 
1=1 aZl 

It is easily verified that the eigenfuctions in this case are functions of the 
type flU), which are homogeneous polynomials of the order k i.e. in (JJ) the 
sum is over k\ + ... + k n = k and the corresponding eigenvalues are Ek = k. 

III. The Transformation to a Canonical Form 

In the Bargmann-Fock representation of a Hilbert space quantum mechanical 
equations for bosons interacting with a manifold of fermion states are rep- 
resented by a system of linear differential equations in the complex domain 
H f|. The physical solutions ([I]) of such a system must belong to T n i.e. 
be entire and obey condition ([|). In practice the solution of this equation is 
complicated due to the occurrence of finite regular singularities. 

In this section we describe a transformation - due to Birkhoff - which 
allows, for a system of linear equations of the first order in one variable, 
to find a canonical form. To be in the canonical form the system must 
be transformed in such a way that all its finite singularities reduce to only 
one singularity at zero, while at the same time preserving the order of the 
singularity at infinity. Consequently, the transformed system is more likely 
to be exactly solvable. 

The system of m linear differential equations of the first order has a 
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general form 

df r m 

-J- = z2Prs( z )fs, r = l,...,m (9) 

az s=1 

and we assume that p rs {z) are analytic functions of a complex variable apart 
from a finite number of regular singularities (even at infinity). Practically it 
means, that outside the circle \z\ = R, which includes all the finite singular 
points, the coefficients may be expanded in a Laurent series 

Prs(z) = £ p^Z k , p« G C, (10) 

k=— oo 

where q > — 1 and q + 1 is termed the rank of the singular point at infinity. 
In general the system @ can have m independent sets of solutions fP(z), 
where s = 1, m denotes different solutions within a set t = 1, m. Now 
we assume a linear transformation of the form 

m 

f r ( Z )=Y / ars(z)F s (z) (11) 
s=l 

where the coefficients a rs (z) are analytic at infinity and reduce at infinity to 
a unit matrix 

oo (k) 

a r s{z) = Y.^ti <&eC, a^ = S rs . (12) 

k=0 Z 

In a sense this {a rs (z)} matrix could be said to contain all the finite singu- 
larities of the initial system. Under this transformation the original system 
d^) turns into a system of a slightly different form 

dFJz) 



z- 



dz ^ , 



J2Prs(z)F s (z), r = l,...,m. (13) 



The coefficients of the transformed system are given by the equation 

{Prs(z)} = z ({a; k \ Z )} {p kl (z)} {a ls (z)} - {a^z)} , (14) 

where {a'^z)} is the matrix inverse to the matrix {a r k(z)} and {a r k(z)}{pk s (z)} 
denotes the matrix multiplication. Now we are in the position to formulate 
the Birkhoff theorem - the crucial result in our analysis. 
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The Birkhoff theorem states that for every system of the type ^) there 
exists a transformation matrix ftilty such that the coefficients P rs {z) of the 
transformed system (|73jj are polynomials of a degree not exceeding q + 1 J7| . 

There are several properties of the above transformation which should be 
stated here. 

(i) All finite singular points of the initial system (||]) coalesce to only one 
singularity at zero (because P rs {z) are polynomials). This is the most signif- 
icant property of the transformation and because of it the system (|H|) may 
be termed the canonical form of the system @ . 

(ii) The ranks of the singular points at infinity for both systems (|9]) and 
are equal. 

(iii) If a given set of solutions f^'(z), s = 1, m, belongs to the Bargmann- 
Fock space T\, then the corresponding transformed solutions F^'(z), which 
are found due to the Birkhoff theorem, also belong to this space. This can 
easily by shown by adapting the treatment by Birkhoff to the case of entire 
functions. 

The final property allows to reject all the solutions of the transformed sys- 
tem which do not belong to the Bargmann-Fock space as being non physical. 
However the inverse of this property is not automatically true, but requires a 
proper choice of the transformation matrix. We will come back to this point 
when discussing the actual solutions of the system under investigation. 

A useful test to check if the solutions of the transformed system can be 
entire is given by the indicial equation of the transformed system || 

det{c rs — pS rs } = (15) 

where c rs = P rs (0), r, s = 1, m. The solutions of (|13|) depend on the roots 
Pi,...,p m of the indicial equation. There are several rules connected with 
this equation which indicate the possibility of existence of entire solutions 
and their degeneracy. 

Dq. If none of the roots is a non- negative integer, the equation (|T3| ) has no 

entire solutions (because they are not analytic in the origin). 

D\. If one of the roots pt is a non- negative integer and the remaining roots 
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are either non- integers or are equal to p t , then there exists exactly one, up 
to linear dependence, set of analytic solutions of and it is of the form 



F^(z) = z pt u s {z), s=l,...,m (16) 

where the u s (z) are analytic and u s (0) ^ 0. Functions (|I|) are entire provided 
their radii of convergence are infinite. 

D2- In the remaining cases if two or more roots of the indicial equation 
are integers (at least one of them is non- negative) , there exists at least one 
analytic set of solutions of the form ([IB]), where p t is the maximal integral root 
of (|15|). The remaining solutions corresponding to other integral roots are 
usually singular in the origin, but in exceptional cases may also be analytic. 

It should be pointed out that the Birkhoff theorem is an existence theo- 
rem, which as such does not provide the actual form of the canonical equation 
nor the transformation matrix. In practice at least the canonical equation 
can usually be found relatively easy by the use of (0), keeping in mind that 
the P rs (z) are polynomials of a given degree. To this aim we invert the equa- 
tion ([H]) and express it in terms of the relevant expansion coefficients. This 
yields, for every integer / > 0, a system of m 2 equations 

e (Kr >} {p£ +1 - i) } - { P r } } Kr } }) = (i-*-i) r-- i} } , (17) 

where expansion coefficients for {p r k(z)}, {a r k(z)} and {P rs (z)} are defined 
by (|TTj), (|TJ) and P rs (z) = E/£o -Prs^ respectively. We assume also that 
{a^l} and {P^ k } with negative indices i are zero matrices. 

This formula can now be used to find the coefficients P rs (z) of the trans- 
formed system. In this case we only need the equations corresponding to 
I = 0, q + 1. Here the trivial case with I = immediately yields: 

{^ +1) } = {^}- (18) 

For higher I, 1 < I < q + 1 the resulting expansion coefficients in the trans- 
formed system may also depend on the aj^ , ...,aj^, +1 ^ coefficients in the ex- 
pansion of the transformation matrix. These coefficients thus may enter the 
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transformed system as extra degrees of freedom, which we will denote as the 
parameters of the transformed equation. The canonical transformation can 
only be defined up to these parameters. However in the context of a physical 
model, their values will be constrained by the requirement that the solution 
of the initial system belong to the Bargmann-Fock space. 

The remaining equations in the system ([H]) i.e. the formulas correspond- 
ing to I = q + 2, q + 3, ... form a set of recurrence equations for the {a rs (z)} 
matrix. This system determines {a rs (z)} as a function of the parameters 
of the transformed system. The procedure of determining the transforma- 
tion matrix for a given set of parameters is in general infinite and it may be 
very difficult to find the parameters that lead to solutions in the Bargmann- 
Fock space. Hence it is conceivable that we know the initial (Q) and the 
transformed (|T3|) systems of equations, without being able to solve the trans- 
formation matrix fll2|) . In some cases this still allows us to draw some im- 
portant conclusions about features of physical interest, such as degeneracies 
or symmetries. 

A special case arises if we assume that the expansion flT2| ) of {a rs (z)} in 
negative powers of z is finite. In this case confinement to the Bargmann- 
Fock space can indeed easily be guaranteed. If the highest order in the 
denominators of (|T2|) is not greater than p t from flI5|), which is the lowest 
power in the expansion of the F®(z), then the corresponding solution of the 
initial system is automatically analytic in the origin. In this case the system 
fll7l) also remains finite and can be solved. This procedure precisely leads to 
the isolated exact solutions of the initial system. 

IV. The Solution of a Two Level System Cou- 
pled to a Single Quantized Mode 

A. The canonical form 

In this section we derive the canonical form of the dynamical equations for 
a two level system coupled to a single quantized mode. The Hamiltonian of 
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such a system, sometimes called Rabi Hamiltonian [|J, is of the form 



H = ua + a + fia 3 + X(a + + a )(a + + a), (19) 

where a + and a are boson field (0) creation and annihilation operators, a ± = 
± «cr 2 ) and a±, o-i-, 03 are Pauli matrices. The parameter u is the boson 
field frequency, 2/i is the atomic level separation and A is the atom - boson 
field coupling constant. We choose the energy unit in such a way that uj — 1 
and we assume that A and \i are not vanishing simultaneously. 

The first step to solve the system fll9f) is to make a unitary transformation 
which replaces operators U\ — > a 3 , a 3 — > and cr 2 — > — cr 2 . Then we write 
the stationary Schrodinger equation for the two-component wave function in 

the position variable ( J • I n ^ ne second step, by replacing a + — > 2 and 

a — > j-, we perform a transition to a Bargmann-Fock space. In this space the 
Schrodinger equation is equivalent to a system of two first order differential 
equations for the Bargmann-Fock space functions fi(z), / 2 (^) € T\ 

d E-Xz ix 

(1Z Z ~T A Z ~v A 

^Hz) = -^- fl ( z ) + ^±^f 2 ( z ), (20) 
az z — A z — X 

where E is an eigenenergy of 7i. Note, that in the Schrodinger representation 
of creation and annihilation operators we have a + — > 4j (^ — ip^), a — > 
-^(^ + «P$), where £ and are conjugate position and momentum. The 
system corresponding to (|20|) in this representation consist of two second 
order differential equations in a real variable £. 

The present, Bargmann-Fock space formulation of the problem has been 
investigated earlier and approximate solutions have been found j|, [§. In 
addition Kus has derived some isolated exact solutions, corresponding to 
degenerate levels 0. 

For our purposes we point out three properties of the system. Firstly, 
it is of the form @ with two finite singularities in z = A, —A. Secondly, 
expanding coefficients of fl20D in a Laurent series fllOD we find that q = and 
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therefore the singular point at infinity is of the first rank. Finally, we note 
the following symmetry: if ^ ^f^(z) J * S a so ^ u ^^ on °^ then ^ ^ 
is also a solution, corresponding to the same energy value. 

The Birkhoff theorem is found to apply to the system in (p0|) . Hence one 
can claim the existence of a canonical form. With q = 0, the coefficients of 
this form will be polynomials of the first degree! They can be obtained from 
equations (|P7D and fllTf ), as explained in the previous section. The linear 
terms of P rs (z) are inferred at once from (|T8|) i.e. P}l'(z) = (— l) r XS rs . To 
calculate the remaining four zeroth order terms of P rs {z) we use equation 
(|T7|) with / = 1. The canonical form of the system (EUf) is therefore 



zj-/i{z) = (E-\z + \ 2 )F 1 (z) + (- f ,-2\a ( 3)F 2 (z) 

z Tz F2(z) = (-^ + 2Xa 2i) F ^ z ) + { E + Xz + x2 ) F ^ z )i ( 21 ) 

where Fi(z) and F 2 (z) are linearly transformed f\(z) and /a (2) (O)' an d 
, are parameters belonging to the transformation matrix (|T^). 
The main feature of the canonical system is that it has only one singularity 

at z = 0, and because of that is exactly solvable. Prior to solving it however 

' AC*) ' 



we exploit the symmetry of solutions I ^ M of the initial system (£(]) • To 



preserve this symmetry in the transformed pair, i.e. if ( ^ ) is a solution 



F 2 (z) 



F 2 (-z) 



of (|2l|) then so is I ^ , I I , we must impose the following symmetry of 



the transformation matrix 

aij(z) = a [i+1][j+1] (-z), i, j = 1,2 (22) 

where [• + •] denotes an addition modulo 2. 

B. Quantization conditions and canonical solutions 



Physical constrains require the solutions of (|20|) to belong to the Bargmann- 
Fock space, i.e. to converge in the entire plane. This leads to quantization 
conditions, as we will show in this section. 
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The canonical system contains four parameters A, /i, E and . In the 
usual perception of the problem A and p are external quantities, describing 
the physics of the system, while E and are essentially free parameters, 
to be determined by the quantum conditions. The indicial equation flTSP 
corresponding to the system (pi) is seen to link i? and : 

p = £ + A 2 ± A, (23) 

where A = p+2Xa^ . For the system (^) to have solutions in the Bargmann- 
Fock space at least one value of p must be a non-negative integer. The value 
of A thus determine the whole energy spectrum but in turn it is controlled by 
the aw parameter, from which the transformation matrix can be generated 
via ([T7|). For the time being we will treat it as a free parameter, delineating 
classes of solutions in the energy spectrum. As we will see later its value will 
be fixed by the requirement that the solutions of the original system (|20|) are 
in the Bargmann-Fock space. 

The general solution of ( 2~ID can be obtained by transforming it into a 
second order form 

, 2 



z 2 F?(z)+z 



2 [E + X' 



F[{z)- 



E + k 



A z + Xz 



AV 



F x {z) = 0. 
(24) 

This equation can be reduced to a confluent hypergeometric (or Kummer) 
equation. Its general solution is a combination of two functions 



F 1 {z) = d exp(Az) 1 F 1 (l + A,l + 2A;-2Xz) z 
C 2 exp(A^) iFi(l-A, 1-2 A; -2Xz) z 



E+X' 2 +A 



+ 



E+\ 2 -A 



(25) 



with arbitrary C\ and C 2 . The function \F\{a, c;z) is called confluent series 
or Kummer function and is defined for all complex a, z and c ^ —n, n = 
0,1,2,... 

„, , a z a(a+l) z 2 a(a+l)(a+2) z 3 , . 



cl! c(c+l) 2! c(c+l)(c+2) 3! 
The Kummer function is entire. The solutions for F 2 (z) are of the same 



general form (25|), with however z replaced by —z. 
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The asymptotic behavior of the solution fl2"5|) for \z\ — > oo is restricted by 
the function z ai exp(Az) + z a2 exp(— Xz), where ai, a 2 are real numbers flQ| 
and therefore the condition (PS) is always obeyed. Consequently the solution 
(P5|) belongs to the Bargmann-Fock space JT X provided at least one of the 
roots of the indicial equation (|23|) is a non-negative integer. 

A particularly simple case occurs if A = because then the transformed 
system (^Tj) coincides with the initial one (|2~0D . The transformation matrix 
reduces in this case to the unit matrix. 

Let us now discuss the solutions fl25|) and their degeneracy for different 
values of E + A 2 . 

(i) If E + X 2 is not an integer nor a half- integer then, to get physical solu- 
tions, we should take A such that one of the numbers p\ = E + X 2 + A and 
p 2 = E + A 2 — A is a non-negative integer. In this case, according to the 
general rule Di, one of the two functions composing Fi(z) can be entire. 
The corresponding solution for F 2 (z) is also one- dimensional and is given by 
F 2 (z) = (-l)^+ t F 1 (-z). 

(ii) If E + X 2 is a half-integer then we take A also half-integral. Despite the 
fact that both E + X 2 + A and E + X 2 — A are now integers (rule D 2 ), The 
solutions F\(z) and F 2 (z) are still one- dimensional because in this case one 
of the numbers 1+2A or 1—2A is a non-positive integer and the corresponding 
Kummer function is not defined. 

(iii) If .E + A 2 is an integer we can take the simple choice A = 0. Then, for 
E+X 2 > each of the solutions is entire and forms a one-dimensional space 
(rule .Dx), but because the system (^l|) is diagonal we can always take its 

two linearly independent solutions of the form ^ F^^(z) J anC ^ 
Consequently in this case the solutions are degenerate. This class of solutions 
will comprise the exact solutions found by Kus. 

(iv) If E + X 2 is an integer and E + X 2 > 0, we can also take A integral, 
< \A\ < E+X 2 . In this case, although the system (|2l|) is no more diagonal, 
solutions Fi(z) and F 2 (z) become two-dimensional and lead to degenerate so- 
lutions. This represents the rarest case of D 2 when two roots of the indicial 
equation are integral and both corresponding solutions are analytical. 
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To summarize we conclude that the only case where solutions are degen- 
erate can occur when E + A 2 is a non-negative integer. 



C. The isolated exact solutions 

In the previous section we have shown that the canonical form of the Rabi 
Hamiltonian can be solved within Bargmann-Fock conditions. This results in 
a coupling between A and E parameters, which is interesting in its own right, 
but does not lead to quantized energies. As we have already mentioned, the 
true quantization condition stems from the requirement, that the solutions 
of the original system belong to the Bargmann-Fock space. This implies that 
the transformation of the canonical solutions must act within the Bargmann- 
Fock space. In this way we fix A and hence E. 

In general this procedure is nontrivial since the transformation matrix is 
generated by an infinite system of equations ([H]), yielding an infinite series 
of coefficients In this section we will not be concerned with the general 
case but only study the exactly solvable class of solutions, which corresponds 
to A = 0. We assume that the transformation is nontrivial i.e. A 7^ 0. 

Let us first consider the simplest possibility \i — 0. With A = 0, and 
hence = 0, the indicial equation yields that p = E + A 2 must be a non- 
negative integer. It is easy to find that the corresponding transformation 
matrix is diagonal and reads 

/(l + Af +A2 \ 

whereas the solutions are generated by 

F 1 (z) = exp(-Xz) z E+x2 

F 2 (z) = exp(Xz) z E+x2 . (28) 



Note that in this case the expansion ([12] ) in negative powers of both diagonal 
terms terminate. The last non zero coefficients are ai E+x2 \ On the other 



hand the lowest power of z in the expansion of the transformed solution (25) 
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is also E + A 2 . It is then easily verified that the solution of the initial system, 
which follows from (|Tl"D, is of form ([[]) and thus belongs to the Bargmann-Fock 



space. The energy spectrum is obtained directly from the indicial equation 
fl2"3j ) and reads E p = p — A 2 , where p = 0, 1, 2, ... . 

In the line of this example we now turn to the more general case /x ^ 
but keeping the a^J matrix finite. As a first example we assume that the 
expansion (|12"D of the transformation matrix terminates, and a^J are the only 
non-zero coefficients. In other words afj = a® = ... = 0. In this case the 
system (|T7|) is also finite because starting from I — 3 all the higher order 
equations are equivalent to equations corresponding to I = 2. Taking into 
account symmetry properties (^2|) of the matrix, the system reduces to four 
equations 

flu + 2pa 12 + 2Aa 2 2 - X(E + A 2 ) = 
2jMn + ai2 + 2Aanai2 + A/i = 
(E + \ 2 )a u + fia 12 - \(E + \ 2 ) = 

pan + (E + \ 2 )a 12 + Xp = (29) 

where we denote ~a~ki = for convenience. Let p, be non trivial and < 
fi < (E + A 2 ). Then, after some algebra, one can show that the system (215) 
is equivalent to the transformation matrix 

/ 1 + mt 
K.(*)}= \t z 

V 2Xz 

and two other conditions: E + A 2 = 1 and 4A 2 + p 




element in (|30|) it also follows that A = p + 2Aa^ 2 

Thus the very assumption that the transformation matrix terminates gave 
us the form of the matrix, the values of E + A 2 and A, showing that they 
correspond to a degenerate solution of (|25|). In addition the system has energy 
E = 1 — A 2 only if the atomic level separation 2p and the atom-boson field 
coupling A obey condition 4 A 2 + p 2 = 1. The corresponding eigenfunction 

( i^jzj ) i s ffl ven by (P8|) an d ( j 1 ?\ ) can be found by applying (]TTJ) . One 
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can easily shown that fi(z) and f 2 (z) belong to the Bargmann-Fock space. 
This solution agrees with the first root of the Kus series ||. 

In an analogous way one can show that if the transformation matrix 
terminates from the third order coefficients onwards yj = ... = and 

< fi < (E + A 2 ) then the transformation matrix functions are 



an(z) = a 22 (-z) = 1 + — + r 2 8A 



z z^ 

k2 , ,,2 



a^z) = a 21 (-z) = -— + — (31) 



i.e. again A = and the remaining conditions are E + A 2 = 2 and 32A 4 — 
32A 2 + 12A 2 /i 2 — 5/i 2 + /i 4 + 4 = 0, which corresponds to the second root of the 
Kus series. In this way all solutions of the Rabi Hamiltonian characterized 
by a terminating {a rs (z)} can be generated. They are found to coincide with 
all known exact solutions. 



V. Discussion 



The Rabi Hamiltonian describes the coupling between a two level system and 
a single mode through a linear interaction term Acr 1 (a + + a). This problem 
is relevant in quantum optics but also appears in molecular physics as the 
simplest example of vibronic coupling |i~2f . 

In Bargmann-Fock space it can be represented by a system of two first 
order differential equations, as shown in (|20|) . For \x = this system separates 
into two independent equations 

(z + \)^-h{z) = {E-Xz)h{z) 
az 

(z - X)^- z f 2 (z) = (E + Xz)f 2 (z). (32) 

This situation corresponds to a highly symmetric case with two uncoupled 
harmonic oscillators that are displaced to the left and to the right in coor- 
dinate space [@, and cross at the coordinate origin. From this perspective 
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the introduction of the level separation 2/i may be viewed as a symmetry 
lowering perturbation, which couples the two oscillators. 

If we now compare the original system (KtJ) to the canonical one (PTj), we 
still have essentially a set of two displaced oscillators - but now the coupling 
term no longer corresponds to [/, but to / u + 2Aa^! Hence the canonical trans- 
formation provides a degree of freedom that allows to perform a displacement 
in the space of the coupling parameter itself. We can thus adjust in such 
a way that it compensates for the energy gap, by requiring — "2A > or 
A = 0. System (pTJ) then becomes 



z-^F 1 {z) = (E-\z + X 2 )F 1 (z) 
clz 

z^F 2 {z) = (E + Xz + \ 2 )F 2 (z). (33) 

This system again describes two degenerate harmonic oscillators that are 
displaced in coordinate space and also underwent a translation in z space 
over + or —A. 

We thus have shown that the exact solutions , that are found for A = 0, 
can be mapped onto the energy spectrum of a degenerate harmonic oscillator. 
A similar result can also be obtained for the Juddian exact solutions of the 
E <S> e Jahn- Teller Hamiltonian IfTH . 
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